A Realistic Cosmology Without Cold Dark Matter 
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We propose a new framework unifying cold darlc matter (CDM) and Modified Newtonian Dynam- 
ics (MOND) to solve their respective problems on galactic scales and large scale structure formation. 
In our framework the dark matter clusters on large scales but not on galactic scales. This environ- 
ment dependence of the dark matter behaviors is controlled by a vector field, which also produces 
the MOND effects in galaxies. We find that in this framework only a single mass scale needs to be 
introduced to produce the phenomena of CDM, MOND and also dark energy. 
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Introduction: The observed universe appears not to 
be purely made of standard model particles and Einstein 
gravity. We are yet to identify the physics of the miss- 
ing constituent (s). Data on galactic and larger scales are 
often used to argue for the MOND and the CDM frame- 
works respectively. Linear growth of large scale struc- 
ture in the early universe, such as the Cosmic Microwave 
Background (CMB), favors the CDM idea, and nonlinear 
structures on the scale of galaxy clusters agree very well 
with numerical CDM simulations. On the other hand, 
on smaller scales, Milgrom's MOND formula [H captures 
the tight correlation of the observed baryonic mass dis- 
tribution within a spiral galaxy vs. the observed gravi- 
tational acceleration at each radii of that galaxy. This 
applies to galaxies with a wide range of scales, forma- 
tion histories and environments, from dwarf to elliptical 
galaxies 0, d, 0i S @ • The amazing accuracy of this re- 
lation and the fact that it predicts correct TuUy-Fisher 
relation even for tidal dwarf galaxies [3] motivate the non- 
covariant MOND theory 181 and a class of covariant the- 
ories fd, 13, [nl, [3 [3, [hIiI], to ehminate the need for 
the CDM particles. 

However, both CDM and covariant MOND have their 
own problems. So far the most challenging difficulty for 
covariant MOND theories is in producing early growths of 
large scale structure and fitting the CMB data [lil. [itI [l8l| 
on which CDM works very well. Also, massive neutrinos 
seem to be indispensable even in covariant MOND to ex- 
plain the lensing data in galaxy clusters 3 lOl ■ In com- 
parison, thanks to its simplicity, the CDM framework en- 
joys many tools for sophisticated numerical simulations, 
yet the properties of galaxies in these simulations are not 
in good agreement with observations. The overproduc- 
tion of dark structures in small scales is well-known as the 
substructure problem and the cusp problem. A common 
assumption to solve this problem is that CDM particle is 
ballistic and will alway follows the same geodesic equa- 
tion as a coUisionless star would. 

Indeed, CDM and MOND arc both complementary 
and mutually exclusive: if both exist in galaxies, then ob- 



TABLE I: The behaviours of the Ci terms in different relevant 
regimes (-1- means there is effect and — means no effect). The 
static limit of the terms actually depend on spatial configura- 
tion of the vector field, but this can be consistently regarded 
as of higher order. 
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viously new problems will arise. A natural way out of this 
dilemma is to have a "CDM" which is no longer cold in 
the environments where MOND dominates, e.g., letting 
it develop a nonzero pressure or have a much smaller mass 
(for particles) there. This environment dependence may 
be controlled by a scalar field, but in this case the differ- 
ent dark matter behaviours in different regimes (galactic, 
cluster and cosmological) indicate that the coupling be- 
tween the scalar field and the dark matter must be fine- 
tuned (if it is possible anyway) , because of the dynami- 
cal nature of the scalar field. In order to see why, note 
that the MONDian behaviour is only expected where the 
Newtonian acceleration |V$| is smaller than the MOND 
parameter ao (on galactic scales), but not in the solar 
system (where |V$| 3> ap), nor on cosmological scales in 
most of the cosmic history (where cH » oq). These sug- 
gest that we should use both |V$|/ao and cH/ao as the 
criteria about the environment dependence of the dark 
matter behaviours. A scalar field dark matter faces not 
only the challenge to reproduce MOND when |V$| < oq 
[2ll |. but also to follow both V$ and cH through a correct 
dynamical evolution, because there are no inherent char- 
acteristic quantities which mimic V<I> or cH in normal 
scalar field models. 

In most attempts to construct relativistic MOND, a 
vector field is used, which can easily overcome the sec- 
ond challenge faced by scalar fields. Furthermore, from 
previous studies of time-like unit-norm vector fields (the 
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iEther field JEa [24I), we know that there are four possi- 
ble kinetic terms for Ma, which have different properties 
in different regimes. If we write these kinetic terms as 
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f^s^d^c + CiM'^I&gcd, in which c,'s are dimensionless con- 
stants, then Table U briefly summarizes the behaviours 
of these terms in different regimes. More specifically, in 
the static limit the ci and C4 terms are oc (V$)^ while in 
background the ci,2,3 terms are oc {cHY (see [ll] for an 
earlier discussion about this point). 

These facts suggest that the jEther field Ma should be a 
natural candidate to control the environment dependence 
of the dark matter's behaviour, most probably through a 
coupling to the latter. In this work we shall give a sim- 
ple example to illustrate this principle. The idea is very 
straightforward: let the behaviour of the dark matter de- 
pend on /C/oq, which could be very large on cosmological 
scales in most of the cosmic history (/C ^ (?H^ ^ Oq) and 
in the solar system (/C ~ |V<l?p S> ag), but becomes of 
order unity or less on galactic scales JC ~ |V$p ^ ao- In 
the first case the dark matter has (nearly) zero pressure 
and zero sound speed, and thus is actually cold, while in 
the latter situations it acquires a nonzero pressure and 
nonzero sound speed, and thus no longer clusters; instead 
then, the Ma field wiU produce the MOND effect. 

The Model: A distinction between our model and pre- 

BEE, 



vious ones 
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18| is that we arc not using Ma 
to grow the large scale structure, which proves difficult. 
Instead, we introduce a dark matter which is controlled 
by Ma through a coupling. For illustration purpose, we 
take (/3 as a k-essence field [i^l , though our principle can 
be applied much more generally (see a discussion below). 
We start from the following Lagrangian density 



C 



-g{v,X,JC) + V{v,]C) + \{M''Ma~l) (1) 



in which g, V are arbitrary functions, is a dimensionless 
scalar field, X = i V^jVotp and /C, Ma are defined above; 
A is a Lagrange multiplier ensuring that Ma has unit 
norm. For g(ip, X, IC), we use a generic power-law in X 
with the power dependent on /C: 



g{^,X,IC) = g(X,IC) = 2w ( ^ 



(2) 



in which the normalization Og is always justified because 
(p (hence X) could be rescaled; w = w (^^^ > which has 
the meaning of the dark matter equation of state, is a free 
function introduced to produce the desired effects, and 
/C3 = c^S^S^iX/ aM'^'^ hM'^ only has the C3 term. Here C3 is 



a constant, and we adopt ^ 3+ j , where n = I for 

illustration. We further choose V{(f, IC) = V (/C4) where 
/C4 = C4{w)M''M'>gcd^aM''\'bM'^ has only the C4 term; 
here C4 can be a function of w, but for the simplicity 
of arguments, we shall first treat C4 as a constant first. 
Note that our choices of the C3 and C4 terms above are 
designed so that we have a clean separation between the 
dark matter behaviours in cosmological background (g) 
and in the static limit {V). This is achieved because the 
C3 term only has effects on the former while the C4 term 
only affects the latter (c.f. Table |T|. 

A variation with respect to the metric tensor gives the 
following energy momentum tensor ^ 



S-KGTab 

= gx'^a^p'^bV - gab[g{X,IC) -V{IC)] 

-2V, [{Vk - gK) {jE'J{ab) + ^(aJ% 
+2Vc [{Vk - gK) J'l MdMaMb 

-2CiVK {M'^VdM,) (^^Ve^") MaMb 

+2ciVK {M^^cMa) {M'^'^dMb) , 
where we have defined J° = K°-\^VbM'^ and 
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gK 



dgiX,JC) _ dV{p,)C) 

Vk = 



dK. 



die 



(3) 



(4) 



We could also derive the scalar and vector field equations 
of motion, but they are not needed here. 

The energy momentum tensor for the k-essence field is 
defined as (see below for a discussion) 



gx^a^^bv - gabg(x, JC) 



(5) 



which resembles the energy momentum tensor of a per- 



fect fluid T^^ = (p^ -t- p^)uaUb " p^gab with Uc 



/2X 



and energy density p^p and pressure : 

SnGpy, 
%ttGp^ 

Substituting Eq. into Eqs. ([51 [T]) it is easy to check 
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2Xgx~g{X,lC), 
g{X,lC). 



(6) 
(7) 
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So we see that when /C 3> a§ this behaves as dust while 
when /C <C flg it behaves as radiation. Furthermore, when 
it behaves as dust the sound speed satisfies 



cl = 



gx 



1 



2Xgxx + gx 







(9) 



^ Note that in Eq. JS} the ?C in g{X, K) can be different from that in 
K) {e.g., having different Ci's). Consequently, if for example 
the K, in V{tp, fC) has no C4 term, then C4V1C is zero, but this does 
not necessarily mean that C4gic is also zero and vice versa. 



SO that it has the desired clustering property of CDM. 

Remember that we have chosen the /C in w ( | as to 
have only a C3 term. It requires no fine tuning to choose 
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C3 SO that all through the cosmic history /C ^ Sc^H'^c'^ 
Oq, since cHq ^ 6ao and H > Hq at earlier times, as the 
result of which a choice of C3 ^ 0(10° — 10^) guarantees 
-7- ~ > IO8C3 1 (if we use n > 1 in w then C3 
could even be set to 1). This indicates that in our model 
the scalar field does behave like CDM in the background 
expansion and the large scale structure formation, where 
there is significant Hubble expansion. 

At this stage one may be worried about the other terms 
in the energy momentum tensor in Eq. ([3|) : are they large 
enough to spoil the good CDM behaviours we have ob- 
tained so far? We discuss how V{IC) and Vjc are negligible 
separately below. For the gjc terms, from Eq. ([2|) we get 



9K = 



2w 



To see that <C 1, note that = wp^ • 
K = 2,cz{cHf > al and 8TrGp^/3{cHf 



2al 



K. 



AirGp^ 



(10) 

Ptpof^/lC when 
■ 0(1), so 

0{l). (11) 



Meanwhile, the current fractional energy density of dark 
matter is 0.2, which means that SirGpipo ~ 0.6(ci?o)^ ^ 



20a§, so we have SnGp^ = 2 



X 



al - 20a^(l + zf 



in which z is the redshift, or j = B{1 + z)'^ with 
10. As a result log ^ = :|^[log B + 3 log(l + z)] - 
0(w) <C 1 today; log (^^^ /w increases with redshift log- 
arithmically at high rcdshifts, e.g., log (^^^ /w ^ 0(10^) 

at z ^ 10^". But w cx (l-|-z)~^ and (l + z)^^ in the mat- 
ter and radiation dominated eras so that indeed both w 
and log (^^) decrease quickly with redshift. This above 

analysis shows that gic ^ O{10''^ — 10^)w ^ 1 in all the 
cosmological epochs of interests, which is easy to under- 
stand because g ^ is very small while ^ is very large 

(this order-of- magnitude estimate holds for general ns). 
The smallness of g/^ strongly suppresses the effects of the 
iEther terms in Eq. ([3]), making them negligible. In fact, 
the in Eq. ([5]) is not conserved, but the smallness of 
gjc implies that the energy exchange between ip and Ma is 

0{w) < 1, Eqs. (35, 



just negligible. Then, asH,f,2£ 



36) of [2JI show that the perturbation growth also mim- 
ics that of CDM for reasonable parameters. Numerical 
results will be reported in a forthcoming paper. 

We next consider the cluster scales, where the observa- 
tions are not compatible with MOND alone but necessar- 
ily incurs a certain amount of dark matter. An example 
is the bullet cluster, in which the offset between the gas 
and dark matter distributions is hard to be explained by 
MOND. These scales generally have not decoupled from 
the background expansion, where according to our model 
the dark matter is still cold. Thus this model has the po- 
tential to explain the observations on cluster scales. 



On galactic scales, the spacetime is more or less static, 
which means that is small enough to make w = ^ 

i, so that the pressure support is strong enough to pre- 
vent any further collapse of dark matter. This eliminates 
the CDM in galaxy systems as we expected, since oth- 
erwise CDM and MOND will coexist, spoiling MOND's 
good fit with data. Now, with the scalar field dark matter 
not clustering and the gic (with only C3 term) having no 
effect in static weak field, it is the role of the V{IC) (only 
C4 term) to produce the MOND effect. To do this, let us 
use the metric ds^ = (1 + 2e$)dt^ - (1 - 2e^)dx''dxi and 
write jE" = iJq + which ae° is the perturbation of 

iE" and e is a small positive quantity. Then up to first or- 
der in e, it is easy to derive that Gqo = — 2<I>'^ — 2didi^ 
where we have used the fact that $ = 5" thanks to the 
absence of anisotropic stresses. For the energy density of 
the fields [c.f. the right hand side of Eq. we already 
know that the first line as well as all gjc terms have negli- 
gible or zero effects, also it is easy to show that up to first 
order in e the last three lines all vanish, while the second 
line reduces to — 2Vi (c4VK;<i>'*) . Defining p = 1 — C4V1C 
the Poisson equation now reads (actually there is also a 
V(JC) on the right hand side of Eq. (fT2l) . but this is like 
a cosmological constant and will not cluster) 



2d,[p{x)d^<S>] = SnGpb 



(12) 



where ph is the local baryon energy density and the argu- 
ment of p{x) is a; = 



= 7 — — where we have 

used /C — — C4|V$p up to O(e^) in the static limit. We 
could choose the form of y(/C) or /i as in 11 1 



1 - p{x) 



V{IC) 



1/(0) (l + f) (1 + f 



.(13) 



where 1^(0) = {-Ci)-^{3aof ■ Clearly the MOND limit 
p{x) —^ X is recovered when x if we choose C4 = —1. 
In general V{JC) serves as a non-uniform dark energy po- 
tential |Tlj , whose local minimum can recover the MOND 
equation, and the background value behaves as a cosmo- 
logical constant far away from galaxies V^(0) = (3ao)^. 
Such a V{0) with C4 = — 1 is however not enough to ac- 
count for the dark energy with SttGpde ~ 81ag, and we 
will come back to this point later. 

In the solar system, again the k-essence field (p does not 
cluster, and the C3 term in gic has no effect. But here the 
MOND effect and cosmological constant effect are both 
suppressed because from Eq. (fT3|) V{IC) 54:X~^al <C 
ag, and p{x) — > 1 — (1 -|-a;/3)^'' ^ 1 in the strong gravity 
regime in which x > O (lO®) ^ 1. In fact the Newtonian 



gravity and the PPN limits are recovered [11 1 . 

Discussion: We want to point out that the model de- 
scribed above is only a very simple one for the Lagrangian 
Eq. 11]). One can also, for example, use the oscillation of 
a canonical scalar field around its potential minimum to 
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provide the dark matter, with the steepness of the po- 
tential depending on /C/ag. In a more phenomenological 
way, we could simply postulate a coupling between dark 
matter particles and the vector field (like the coupling 
with a scalar field) as a result of which the dark matter 
particle mass depends on -Sr. Furthermore, it is also in- 

teresting to see if the parameter Oq is indeed determined 
dynamically. These possibilities will be considered in de- 
tails in forthcoming papers. 

The interesting fact ao ~ cHq suggests that there may 
be some fundamental relations between MOND and dark 
energy. In fact, there are many possible ways by which 
our model can be generalized to include dark energy as 
well. One way is that at late times when H ~ Hq the dark 
matter decays into dark energy {e.g.^ its equation of state 
w becomes —1). The idea here is to use the quantity gq 
to determine both the transitions from CDM to MOND 
and from CDM to dark energy. A more straightforward 
method is to have a cosmological constant in V{IC): as 
is shown above, the MOND effect only depends on Vjc 
but not V^(0), and we can use the dark energy density 
to fix V{0) so that the combination of dark energy and 
MOND completely determines V{IC). Another interest- 
ing possibility is to note that in Eq. (fT3|) the MOND effect 
requires C4 = — 1 while dark energy requires C4 ~ — 
This can be easily achieved, again using our principle of 
environment dependence: let C4 depend on /C3, for exam- 
ple with C4 = — (3 — 6w)~^ . In this case V{JC) acts as 
an environment dependent cosmological constant, which 
accounts for the cosmic acceleration in background cos- 
mology {w 0) and approaches zero in the solar system 
{w 1/3). Note that dark energy and MOND are uni- 
fied with a single V{IC) in the latter two possibilities, and 
we thus have a full Lagrangian as 



C = —2wan 



X 



V{JC) 



(14) 



and 



in which V{K.) = (9ao)2 (1 - 2w)' (1 + + f )2 

X — y/JCi/ttQ for the third possibility. Interestingly, 
flo ~ cHq is a single mass scale introduced for this model 
to relate CDM, MOND and dark energy together. All the 
other parameters (qs, n) are dimensionless and ~ 0{1), 
and there are no fine-tunings of them: the huge difference 
between the dark matter density at earlier times and 

the scale comes as a generic result of the dynamical 
evolution of the vector field. In this sense the vector field 
acts as a leverage, making the tiny mass scale qq capable 
to characterize the large energy density of dark matter. 
Meanwhile, this could also shed further light on the dark 
energy coincidence problem, since the dark energy domi- 
nance begins at the time when galaxies have formed (and 
we observers come into existence) , both characterized by 
our fundamental mass scale oq. 

Summary: In this work we have tried to tackle the 
problem of how to unify CDM and MOND in a consis- 



tent way. The idea is to give the dark matter an en- 
vironment dependence, making it behave like CDM on 
large scales, while reproducing the MOND (Newtonian 
dynamics) in the static and weak (strong) field limits 
respectively. Although the idea of an environment de- 
pendence is not new, it is novel to use the vector (jEq) 
field as the switch. We show how the particular prop- 
erties of the vector field make it very effective for this 
purpose. Our model provides a general framework which 
can potentially solve the problems of CDM on galactic 
scales and of MOND on larger scales. It could also be 
generalized to include dark energy in a way such that all 
the phenomena of CDM, MOND and dark energy are re- 
lated to one parameter gq, which is the single mass scale 
introduced in our model. Both fields in Eq. (|14p . likely 
effective, should provide insights to people seeking the 
fundamental fields in particle physics theories. 
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